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I Abstract 
(N ' 

I V-type asteroids are bodies whose surfaces are constituted of basalt. In the Main 

Asteroid Belt, most of these asteroids are assumed to come from the basaltic crust of 
asteroid (4) Vesta. This idea is mainly supported by (i) the fact that almost all the 
known V-type asteroids are in the same region of the belt as (4) Vesta, i.e. the inner 
^ ■ belt (semi-major axis 2.1 < a < 2.5 AU), (ii) the existence of a dynamical asteroid 

, family associated to (4) Vesta, and (iii) the observational evidence of at least one large 

^ ' craterization event on Vesta's surface. One V-type asteroid that is difficult to fit in this 

scenario is (1459) Magnya, located in the outer asteroid belt, i.e. too far away from (4) 
^ ' Vesta as to have a real possibility of coming from it. 

' The recent discovery of the first V-type asteroid in the middle belt (2.5 < a < 2.8 

AU), (21238) 1995WV7 (Binzel et al. 2006, DPS 38, #71.06; Hammergren et al. 2006, 
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[astro-ph/06094 20), located at ~ 2.54 AU, raises the question of whether it came from 
(4) Vesta or not. In this paper, we present spectroscopic observations indicating the 
existence of another V-type asteroid at ~ 2.53 AU, (40521) 1999RL95, and we investigate 
the possibihty that these two asteroids evolved from the Vesta family to their present 
orbits by drifting in semi-major axis due to the Yarkovsky effect. The main problem 
with this scenario is that the asteroids need to cross the 3/1 mean motion resonance 
with Jupiter, which is highly unstable. 

Combining numerical simulations of the orbital evolution, that include the Yarkovsky 
effect, with Monte Carlo models, we compute the probability of an asteroid of given di- 
ameter D to evolve from the Vesta family and to cross over the 3/1 resonance, reaching 
a stable orbit in the middle belt. Our results indicate that an asteroid like (21238) 
1995WV7 has a low probability (less than 10%) of having evolved through this mecha- 
nism due to its large size (I) ~ 5 km), because the Yarkovsky effect is less efficient for 
larger arteroids. However, the mechanism might explain the orbit of smaller bodies like 
(40521) 1999RL95 (Z) ~ 3 km), provided that we assume that the Vesta family formed 
> 3.5 Gy ago. We estimate the debiased population of V-type asteroids that might exist 
in the same region as (21238) and (40521) (2.5 < a < 2.62 AU) and conclude that about 
10% or more of the V-type bodies with D > 1 km may come from the Vesta family by 
crossing over the 3/1 resonance. The remaining 90% must have a different origin. 
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1. Introduction 



The recent discovery of (21238) 1995WV7 iBinzel et alJbnnd i lMasi et al.ll2nnd : lHammergren et al 



20061 ). a basaltic asteroid in the middle Main Belt (2.5 < a < 2.82 AU), raised new questions about 
the origin of basaltic material in the asteroid belt. 

Basaltic asteroids show a spectrum characterized by the presence o f a deep absorption band cen 



tered at ~ 0.9 ^m, and are classified as V-type in the usual taxonomies (jTholen 



1989 



Bus and Binzel 



2002al ). Most of the known V-type asteroids are fragments from the crust of asteroid (4) Vesta. 
This is supported b y the fact that (4) Vesta is the only large asteroid showing a basaltic crust 
|McCord et al.lll970l ). and almost all V-type asteroids are found in the same region of the Main Belt 
as Vesta, i.e. the inner belt (o < 2.5 AU). Moreover, (4) V esta has associated a dyn amical asteroid 
family, the Vesta family, whose members are also V-tvp e (IMothe-Diniz et al. l2005ll and orig inated 
from the excavation of a large crater on Vesta's surface (IThomas et al.lll997l : lAsphaug)ll997l ). 



The discovery of (1459) Magnya |Lazzaro et al.l 120001 1. a V-type asteroid in the outer belt 
(a > 2.82 AU), provided the first evidence for another possible source of basaltic asteroids in the 
Main Belt. (1459) Magnya is too far away from the Vesta region as to have any chance of being 
a fragment from Vesta's crust. No dynamical mechanism is known to be able to bring an asteroid 
from the Vesta family to the Magnya region. Moreover, Magnya is too big (diameter D ^ 1 7 km ) 
as to fit within the size distribution of the Vesta family (see Sect. H]). iMichtchenko et al.l (120021 ) 
suggested that Magnya is the fragment from a differentiated parent body that broke up in the outer 
belt, but up to now no oth er V-type asteroids hav e been confirmed in th e same region of the belt 
to support this hypothesis (IDuffard and Roig 120071 : iMoskovitz et al.ll2007l ). 



The case of (21238) 1995WV7, the first V-type asteroid discovered in the middle belt, has 
some similarities with the case of Magnya, but also shows some differences. (21238) 1995WV7 is 
also far away from Vesta, and ejection velocities larger than ~ 2 km/s would be necessary to directly 
transport it form Vesta's surface to its present orbit. These ejection velocities can not be produce d 
in typical craterization events similar to the one that originated the Vesta family (|Asphaug)ll997l ). 
On the other hand, (21238) 1995WV7 is close to the outer border of the 3/1 mean motion resonance 
with Jupiter (hereafter J3/1 MMR), centered at 2.5 AU, and its size (Z) ~ 5 km) fits within the 
size distribution of the Vesta family. Since the inner border of the J3/1 MMR is very close to the 
outer edge of the Vesta family, the possibility of (21238) 1995WV7 being a former member of this 
family that reached its present or bit after crossin g over the resonance cannot be totally ruled out. 
The J3/1 MMR is highly chaotic |Wisdomlll982l ) and is considered a difficult-to-cross barrier, but 
up to now no detailed studies have been made to confirm this. 



A different scenario is proposed by ICarruba et al.l (120071 ). in which the source for (21238) 
1995WV7 and other, yet undiscovered, basaltic asteroids in the middle belt could be asteroid (15) 
Eunomia located at a ~ 2.65 AU. This asteroid appears to b e a partially differen tiated parent 
body showing a basaltic-like composition in part of its surface (INathues et al.ll2005l ). Carruba et 
al. suggest that several collisions made (15) Eunomia lose its basaltic crust almost completely, and 
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the subsequent fragments were significantly dispersed in the middle belt over the age of the Solar 
System. 

The aim of this work is to study the possible origin of (21238) 1995WV7 and other V-type 
asteroids in the middle asteroid belt. In particular, we analyze the possibility tha t asteroids from the 



Vesta family increase their orbital semi-major axis due to the Yarkovsky effect (IVokrouhlicky et al 



2OOOI ) and cross over the J3/1 MMR reaching stable orbits in the middle belt. In Section O we 
introduce the population of V-type asteroids observed by the Sloan Digital Sky Survey in the 
middle belt, and present our spectroscopic observations with the Gemini South Telescope that allow 
to confirm another V-type asteroid in the same region as (21238). In Section [3j we describe our 
simulations, analyze the resonance crossing mechanism proposed above, and evaluate its efficiency 
to produce V-type asteroids beyond 2.5 AU. In Section [S, we discuss our results compared to the 
debiased distribution of V-type asteroids in the middle belt. Finally, Section [6] is devoted to the 
conclusions. 



2. V-type asteroids in the middle belt 



The existence of V-type asteroids in the middle belt was first predicted by lRoig and Gil-Hutton 



|2006l ) , who analyzed th e colors of the Sloan Digital Sky Survey Moving Objects Catalog (SDSS- 



MOC: llvezic et al.ll200ll ). These authors studied a sample of 13 290 asteroids contained in the 3rd 
release of the SDSS-MOC that show errors smaller than 10% in all the five bands of the SDSS 
photometric system, named u,g,r,i,z, respectively. They found three candidate V-type asteroids 
in the middle belt, which are listed in Tabled! Two of them are located very close to the outer border 
of the J3/1 MMR. The third one is close to the outer border of the J8/3 MMR, ce ntered at 2.7 



AU. T he s pectroscopic confir mation of the basaltic nature of (21238) was reported by iBinzel et al 



(12006!) and 



Masi et al 



(120061). ba sed on spectroscopic observations in the near infrared (NIR), and 



also by iHammergren et al.l (|2006l ) based on visible spectrocopic plus NIR photometry. 



Table 1: Proper semi-major axis Op, proper eccentricity Cp, sin of proper inclination sinip, absolute 
magnitude H and diameter D of predicted V-type asteroids in the middle belt according to the 
colors of the SDSS-MOC. Di ameters have been estimated assuming and albedo of 0.4, similar to 
the albedo of (4) Vesta fO.42. lTedescnlll989l v 



Name 


ap [AU] 




sin/p 


H 


D [km] 


(21238) 1995WV7 


2.54108 


0.1371 


0.1866 


13.04 


5.15 


(40521) 1999RL95 


2.53111 


0.0458 


0.2159 


14.36 


2.80 


(66905) 1999VC160 


2.74627 


0.1457 


0.2291 


14.51 


2.62 



As part of an observationa l campaign to confirm the taxonomy of V-type candidates identified 
bv iRoig and Gil-HuttonI (l2006l l. we obtained spectra of (21238) and (40521) in the visible. The ob- 
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servations were carried out during the nights of 29-30 April, 2006, at the Gemini South Observatory 
(GS), using the Gemini Multi-Object Spectrograph (GMOS). Table [2] provides the observational 
circumstances of the targets. In order to remove the solar signature from the asteroids spectra, we 
also obs erved the stars SA 107-871 {V = 12.4) and SA 110-361 {V = 12.4), taken form the Selected 
Areas of Landolt (|l992l ). that we used as solar analog star^. 



Table 2: Observational circumstances of (21238) and (40521): A is geocentric distance in AU, r is 
heliocentric distance in AU, (j) is solar phase angle, 6 is solar elongation, V is visible magnitude. 
Dates correspond to the starting time of the observations. The observing conditions were better 
than 85% of image quality, 70% of sky transparency (cloud cover), 100% of sky transparency (water 
vapor), 80% of sky background, air mass < 1.5. 



Asteroid 


Date [UT] 


a (J2000) 


5 (J2000) 


A 


r (f) 





V 


(21238) 


2006 Apr 29.4094 


20'^44'"53.0P 


-24°57'12.8" 


2.599 


2.818 20.9° 


91.9° 


18.3 


(40521) 


2006 Apr 30.2811 


14i>ll'"29.86" 


-24°01'44.8" 


1.489 


2.487 4.1° 


169.8° 


18.1 



Tracking of the asteroids at non-sideral rate was not possible because the use of the peripherical 
Wavefront Sensor (WFS) is not recommended due to flexure within GMOS. Instead, we used the 
On-Instrument Wavefront Sensor (OIWFS) tracking at sideral rate, with the slit oriented in the 
direction of the asteroid's proper motion. 

All the observations were performed using the following GMOS configuration: grating R400, 
filter OG515 to avoid second order spectrum contamination longwards of ~ 0.7 /im, slit width 
1.5 arcsec (the maximum allowable with GS-GMOS), central wavelength 0.73 /xm, spectral coverage 
0.522 to 0.938 /im, spatial binning 2, and spectral binning 4. This configuration provides a resolution 
R ~ 3 000 at 0.90 /xm, but R ~ 200 is enough to detect the deep absorption band longwards of 
0, 75 fim typical of V-type spectra. Therefore, it was possible to do a 15:1 rebinning of the asteroids 
spectra to improve the final signal-to-noise (S/N) ratio by a factor of ~ 4. 

Each asteroid was observed six times at six different positions along the slit separated by 10 
arcsec. Each solar analog star was observed three times at three different positions along the slit 
with the same separation. The integration times for the asteroids allowed to attain S/N~ 20 at 
0.90 /im which, after spectral rebinning, provided a minimum S/N ~ 70. Table [3] summarizes the 
asteroids exposure times. The exposure times for the solar analog stars were chosen to reach 50-70% 
of the detector full well (100k electrons) per exposure. Wavelength calibration of the spectra was 
performed using a standard CuAr lamp. 

The GMOS IRAF package was used to perform the standard reduction tasks. Since each image 
was dithered along the spatial direction, we remove the fringes by making a background image 



^These are G2V stars with solar colors in the magnitude range allowed by Gemini telescopes {V > 11). We recall 
that no spectroscopic solar analog stars are available in the literature in this magnitude range. 
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Table 3: Number of exposures, exposure times (individual and total), S/N at 0.90 ;um of the raw 
and rebinned spectra, and solar analog used in the reduction. 



Asteroid 


^exp 


T 

exp 


^total 


S/N (raw) 


S/N (rebinned) 


Solar analog 


(21238) 


6 


200 sec 


1200 sec 


20 


80 


SA 110-361 


(40521) 


6 


500 sec 


3000 sec 


18 


72 


SA 107-871 



resulting from the combination of three successive images centered around the time of the image 
we want to correct. After the fringe correction, the individual frames were coadded to improve the 
S/N and the spectrum was extracted. 

The final spectra of (21238) and (40521) are shown in Fi gure [H We have plot in gray sev- 
eral spectra of know n V-type asteroids taken from the SMASS (IBus and Binzelll2002bl ) and S^OS^ 
(jLazzaro et al.ll2004l ) surveys for comparison. Both spectra are compatible with the V class, showing 
the typical absorption band with a minimum at 0.90 /xm. Our observations with Gemini constitute 
another independent confirmation of the basaltic nature of (21238), and also allow the confirmation 
of (40521) as the second V-type asteroid found in the middle belo 



3. Possible origin of V-type asteroids in the middle belt 



In the following, we analyze the possibility that (21238) and (40521) were fragments of (4) Vesta 
that evolved to their current orbits by the Yarkovsky effect. The main problem with this scenario 
is that these asteroids needed to cross the J3/1 MMR which is notably unst able. To estimate 
the time Atf^rnss that a 5 krn asteroid would require to cross the J3/1 MMR, Masi et al. (|2006l ) 
and Hammergren et al. (|2006l ) divided the resonance width Aa by the drift rate da/dt due to the 
Yarkovsky effect, and concluded that Atcross ^ tmst, where tinst is the instability time scale of the 
J3/ 1 MMR (i.e., the tirn e required for a population of asteroids to be removed from the resonance. 



Gladman et al.l 119971 ). These authors thus inferred that it would be impossible for an asteroid 



see 

like (21238) to cross the J3/1 MMR. However, this argument is only approximative because it does 
not take into account the resonance dynamics. Here we used precise numerical simulations to study 
whether the crossing is actually possible. 

We performed a series of simulations of the evolution of test particles representing selected 
members of the Vesta dynamical family. To define this family, we used the database of 188 207 
asteroid proper elements rele ased by March 200^. We applied the Hierarchical Clustering Method 



(HCM - IZappala et al.lll990l ) and defined the Vesta family at a cutoff of 60 m/s, which is 5 m/s 



^The ASCII files of the spectra are available at |http : / / staff . on.br/froig/vtypes l 

^AstDys 'http: //hamilton.dm.unipi . it/cgi-bin/astdys/astibo' This database is contemporary to the 3rd. 
release of the SDSS-MOC, therefore the two datasets can be directly compared. 
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larger than the predicted quasi-random leve0 for the inner belt. This guaranteed that we defined 
the Vesta family with the largest possible number of members (~ 9 500) that could be detected 
from the given dataset of proper elements. 

Among the members of the Vesta family, we selected 21 asteroids with 2.485 < < 2.490 
AU, which are the closest ones to the J3/1 MMR, and generated 100 clones of each. All 100 clones 
had the same orbital parameters as the original asteroid, but we allowed each to drift in semi- 
major axis at a slightly different rate da/dt > due to the Yarkovsky effect. This implied that 
the clones reached the border of the J3/1 MMR at slightly different phases of the resonant angle 
a = 3Aj — X — w, where Aj,A are the mean longitudes of Jupiter and the asteroid, respectively, 
and Tu is the longitude of perihelion of the asteroid. Therefore, they sampled different resonant 
interaction histories. 

The simulations were performed using a modified version of the SWIFT_MVS integrator 
(jLevison and Duncanlll994l ). The bodies were considered as massless particles subject to perturba- 
tions from all planets except Mercury, and the Yarkovsky effect was introduced in the simulation 
as an additional acceleration term depending on the velocity as: 

d'^f\ _ GM da v 
It ^ 



dt^ ) 



Yarko 



2o2 



where G is the gravitational constant, M the mass of the Sun, a the osculating semi-major axis of 
the orbit, f,v the instantaneous position and velocity of t he body, and da/dt t he required drift rate 
measured in AU/y. According to the analytical theory of IVokrouhlickyI (119991 ). the Yarkovsky drift 
rate approximately scales with diameter as: 



da 
lit 



2.5 X 10~ 



1 km 
D 



cos e 



where D is in km, e is the spin axis obliquity, and the coefficient was obtained assuming physical and 
thermal parameters typical of basalt and albedo 0.4. We chose the thermal parameters to produce 
large but plausible drift rates, because this would favor the J3/1 MMR crossing. Slower drift rates 
would apply if the real thermal parameters have different values. The so-called seasonal Yarkovsky 
effect |Rubincamlll995l ) was not included in our model because it only produces da/dt < and it 
is an order of magnitude smaller than th e diurnal eff ect modeled by Eq. ([J) for km-size aste roids. 
The collisional reorientation of spin axes (|Harri£lll979l ) and the YORP effect (|Rubincamll2000l ) were 
not taken into account either. T he former effect ha s proven not to be so relevant in changing the 
spin axis obliquities as believed (jBottke et al.ll2006l ). while the YORP effect statistically tends to 
align the spin axes around e ~ 0, vr, thus forcing a maximum Yarkovsky drift, on average, for all the 
particles. 

We perfomed two series of simulations. In the first series, we assumed that each clone drifted 
outwards at a rate randomly chosen in the interval (1.0 ± 0.1) x 10^^ AU/My, that according to 



■^The average minimum distance between pairs of neighbor asteroids. 
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Eq. ([T]) corresponds to -D ~ 2.5 km V-type asteroids. These simulations spanned 150 My. In the 
second series, the clones drifted at a rate randomly chosen in the interval (1.0 ±0.1) x 10^'^ AU/My, 
which corresponds to D ~ 250 m asteroids. These simulations spanned 25 My. The time step of 
the integrator was set to be 0.04 y. 

In both simulations, more than ~ 97% of the particles were discarded because they entered 
the J3/1 MMR and chaotically evolved to planet crossing orbits in a few My. These particles were 
eliminated by close encounters with the planets, mainly the Earth and Mars, or by impacting the 
Sun. But a small fraction -less than ~ 3%- of the particles crossed over the J3/1 MMR and ended 
the simulation in stable orbits with a > 2.5 AU. 

Figure [2] shows two examples of test particles that crossed over the J3/1 MMR. The left panels 
correspond to a slowly drifting particle and the right panels correspond to a faster drift. The 
particles enter the J3/1 MMR at the side of lower semi-major axes, perform a few librations around 
the resonant semi-major axis, and exit the resonance at the side of larger semi-major axes. The 
particles remain in the resonance for at most a few 10^ y. Their eccentricities and inclinations are 
not significantly affected by the passage through the resonance. In the examples shown in Fig. [2j 
the eccentricities decay to lower values but this happens after the particles have already crossed 
the J3/1 MME@. In general, the particles that jumped the J3/1 MMR ended the simulations either 
with higher or lower values of the eccentricities and inclinations than their original values. 

In Fig. Owe show the trajectory of the slowly drifting particle of Fig. [2] at the exact moment 
in which it crosses the J3/1 MMR. We can see that the particle spent only ~ 5 000 y inside the 
resonance (top panel), so it virtually "jumps" over the resonance. The numbers in the bottom 
panel provide the temporal sequence of the trajectory. The particle enters the resonance somewhere 
between "1" and "2", then it performs one and a half libration sticking close to the resonance 
separatrix, and finally exits the resonance somewhere between "8" and "9". This trajectory is a 
typical example of the resonance crossing mechanism. We expect that only the orbits that remain 
close to the separatrix would be able to exit the resonance pushed by the Yarkovsky drift, and this 
might depend on the particular phase of the resonant angle, a, that the orbit has when it enters 
the resonance. 

Table 4: Fraction /cross of asteroids that crossed the J3/1 MMR during our simulations for different 
Yarkovsky drifts. 



da/dt xlQ-"^ [AU/My] 


/cross 


1.0 ±0.1 


0.0033 ± 0.0019 


10.0 ±1.0 


0.029 ± 0.005 



^This effect is probably due to the interplay with non linear secular resonances located close to the border of the 
J3/1 MMR, like the g — 2ga + gs resonance, where g is the secular frequency of the asteroid perihelion, and g5,ge are 
the frequencies of the perihelia of Jupiter and Saturn, respectively. 
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From our simulations, it was possible to estimate the fraction of test particles that crossed the 
J3/1 MMR, /cross, and its formal uncertainty, as a function of the Yarkovsky drift rate. The results 
are given in Table [H Assuming that these results reflect a linear dependence between /cross and 
da/dt within the studied range of drift rates, and taking into account Eq. ([T]), we find that 

fcrossiD) ~ 0.0075 ^ cose. (2) 

The application of this formula in the size range 1 km ^ D < 5 km requires only modest extrapola- 
tions of the linear dependence toward smaller values of the drift rate (~ 0.5 x 10~^ AU/My). This 
result demonstrates, for the first time, that the crossing of Vesta family members over the J3/1 
MMR is possible in practice. The crossing probability ranges from 0.15% to 0.25% for asteroids in 
the size range of (21238) and (40521). 



4. The production of V-type asteroids with a > 2.5 AU 

In the preceding section, we have shown that it is possible for a former Vesta family member 
to cross over the J3/1 MMR, reaching a stable orbit in the middle belt. We have also estimated 
the crossing probability as a function of the size (Eq. [2]) to be ~ 0.1-0.8% for km-size asteroids. 
We may now wonder whether this may actually explain the presence of (21238) and (40521) in the 
middle belt. 

The number of Vesta family members, with diameters between D and D + dD, that cross over 
the J3/1 MMR during the age of the Vesta family, Tage, and end in the middle belt is: 

dNa>2AD) = dNo{D) f,ol{D) /reach (I?) /cross (1?) (3) 

where dNQ{D) is the number of Vesta family members that were produced when the family formed. 
The factor fcoi{D) accounts for the fraction of Vesta family members in the interval [D,D + dD] 
that, on one hand, survived and, on the other hand, were created due to collisional evolution in Tage- 
It is worth noting that the product dNQ{D) /coi(-D) is equivalent to the presently observed number 
of family members, dN(D). The factor /reach (-C*) accounts for the fraction of Vesta family members 
that were able to reach the border of the J3/1 MMR in Tage aided by the Yarkovsky effect. Finally, 
the factor /cross(-C) is provided by Eq. ([2]). In order to apply Eq. ([3|), we have to address three 
main issues. 



The first issue concerns the actual age of the Vesta family, Tage- Initial estimates by lMarzari et al 



(|l996l ). based on collisional evolution models aiming to reproduce the presently observed size fre- 
quency distribution of the Vesta family, indicated Tage ~ 1-2 Gy. However, these estimates have 
large uncert ainties due to u n know n initial conditions and to uncertainties in the various collisional 



parameters. ICarruba et al.l (|2005l ). based on the dynami cal evolution of some Vesta family "fugi 



fives", found that Tage ^ 1-2 Gy. On the other hand, Bogard and Garrison (|2003l ). measuring 



isotopic abundances in the Howardite-Eucrite-Diogenite (HED) meteorites -that presumably come 



- 10 - 



from Vesta, (IMigliorini et alJll997l )-. concluded that Vesta's crust suffered several major crateriza- 
tion impacts > 3.5 Gy ago. This could imply Tago ^ 3.5 Gy. In view of the wide range of possible 
ages, we will adopt for our calculations three different values of Tage: 1.5, 2.5 and 3.5 Gy. 

The second issue to address is to estimate the fraction of family members, /reach, that would 
be able to reach the J3/1 MMR in Tage due to the Yarkovsky effect. For this purpose, we used 
a Monte Carlo algorithm. We generated an artificial Vesta family constituted of 10 000 fragments 
with individual ejection velocities, Ve^, attributed by assuming a Maxwellian distribution with mean 



ejection velocity fej2l- We considered only the fragments with 



^'ej - ^esc > and f ej - f < 



where fesc = 314 m/s is the e scape velocitv f rom the surface of (4) Vest80, and Vcnt = 600 m/s 
is a maximum cutoff velocity ([Asphaug) 119971 ) . We further assumed that, at the moment of the 
impact, (4) Vesta had a true anomaly and perihelion argument such that the fin al distribution of 



the f ragments in proper elements space is spread over a wide range of Op, €„, /„ (jMorbidelli et al 



19951 ). Finally, we considered that all the fragments had the same diameter -i.e., the ejection 
velocity does not depend on size-, and attributed to each fragment a random value of cos e between 
—1 and 1. Using Eq. ([T]), we computed for each fragment the total drift over Tage and determined 
the fraction /reach that ended with a > 2.5 AU. 

Figure m shows /reach(-C') as a function of Tage and Vey For very small sizes, /reach weakly 
depends on Vej and tends to 50% for D — > 0. This means that about half of the smallest fragments 
will eventually reach the J3/1 MMR (the other half has cose < so is drifting inwards). On the 
other hand, for very large sizes /reach — > since large bodies are much less affected by the Yarkovsky 
effect. It is worth noting that size range where the behavior of /reach is more critical, i.e. more 
sensitive to the different parameters and especially to the family age, is between 2.0 and 7.0 km, 
which is precisely the size range of the asteroids listed in Table [H 

The third issue we need to address concerns the determination of the size frequency distribution 
(SFD), dN{D), of the Vesta family. The Vesta family, as defined by the HCM (Section [3]), shows 
a very peculiar SFD. Going from larger to smaller sizes, we first find (4) Vesta with D ~ 460 km. 
Then we find four asteroids with 25 < D < 130 km and twelve asteroids with 8 < D < 15 km, which 
taxonomic types, based on spectroscopic observations ( Bus and Binzelll2002al : iLazzaro et al.ll2004l ). 
are different from the V-type. These asteroids are interlopers in the family and we can exclude them. 
Finally, we find a huge amount of members with D < 8 km. This group includes all the vestoids, 
i.e. the known V-type asteroids that are members of the Ve sta family. According to simulations of 
asteroid fragmentation using hydrocodes (jPurda et al.ll2007l ). this kind of SFD, where there is only 
one large asteroid and a huge amount of very small fragments with no bodies in between, would be 
typical of craterization events as the one that presumably formed the Vesta family. 



^ We recall that the specific energy of the impact that generated the family is oc v^j . 



being R the radius of Vesta and p its density ijPetit and Farinellalll993l ) 
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Now, dN{D) is related to the differential SFD, n{D), through 

dN{D) = n{D) clD 

and the cumulative SFD, i.e., the number of family members with diameter > D, is given by 

N{D) = / n{D)dD. (4) 
Jd 

where Dmax is the size of the largest family member (i.e. (4) Vesta). Figure [5] shows the cumulative 
SFD of the Vesta family after removing the known interlopers. For D < 8 km the cumulative SFD 
can be fit to power laws of the form NqD'^ (dashed lines in Fig. E]), where 7 and A^o adopt different 
values at different size ranges. This change in 7 and A'^o is produced by two different effects: 



1. The fact that the sample of known family members is complete only up to a given size. It 
is usually assumed that this completenes s limit correspond s to the size where the cumulative 
SFD shows the first inflexion point (e.g. iTanga et al.lll999l ). which is 5.5 km in Fig. [5l 



2. The natural dynamical/collisional evoluti on of the family rnembe rs, which is known to produce 
a shallow cumulative SFD at small sizes (jMorbidelli et al.ll2003l ). 



Since both these effects are difficult to estimate, we will adopt here two extreme approaches to 
model the cumulative SFD: 



• To use the observed cumulative SFD as shown in Fig. [3 This model partly accounts for 
effect ^2, but has the drawback that it will be biased by the incompleteness of the sample for 
D < 5.5 km. 

• To use a single power law fit between 8 and 5.5 km and to extrapolate it to smaller sizes. 
This model partly accounts for effect #1, but has the drawback that it will significantly 
overestimate the number of small members in the Vesta family, especially in the size range 
2.0-5.5 km that is crucial for our study. With this approach, the differential SFD will also be 
modeled by a single power law of the form no-D^~^. 



Knowing all the quantities involved in Eq. jS]), it is now possible to determine the number of 
Vesta family members with size > D that cross over the J3/1 MMR and reach the middle belt in 



rD' 

J ^max 

Na>2AD)= n{D)f,,,,i,{D)Uoss{D)dD. (5) 

Jd 

This integral can be solved numerically taking into account that -D^ax — ^ km. The results are 
shown in Fig. [6] for Vej = 200 m/s and three different values of Tage- The full lines were obtained 
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using the observed SFD, and the dashed hnes using the single power law model of the SFD. The 
true value of Na>2.^{D), for a given Tage, should be something between the corresponding full and 
dashed lines. From Fig. [6l it is clear that the predicted Na>2.5 for > 5 km is at least two orders of 
magnitude smaller that required to produce a body like (21238). Therefore, we may conclude that, 
even if the crossing over the J3/1 MMR is dynamically possible, it is highly improbable that (21238) 
had reached its present orbit via this mechanism. On the other hand, for D > 3 km the predicted 
-^a>2.5 is compatible with the presence of (40521) in the middle belt, provided that Tage ^ 3.5 Gy. 



5. Predicted vs. observed number of V-type asteroids in the middle belt 

Our results above have been compared to the observed population of V-type asteroids in the 
middle belt, which is most likely incomplete. Here we will estimate the debiased SFD of V-type 
asteroids in the region 2.5 < a < 2.6 AU, where (21238) and (40521) are found, and will compare 
it to our estimates of Na>2.b- 

In terms of absolute magnitude, the unbiased SFD of V-type asteroids in the region of interest, 
n\{H) dH , is related to the observed SFD, ny^{H)dH, through 

n^^" {H)dH = b{H)ny{H)dH (6) 

where b{H) is a bias function that we need to determine. This bias function arises from two main 
effects: (i) the fact that the known population of asteroids is complete only up to a certain size, 
and (ii) the fact that the observed population n'^' is obtained from the SDSS that mapped only 
a fraction of the total known population of asteroids. We may assume that effect (i) affects the 
asteroid populations at both sides of the J3/1 MMR approximately in the same way, so it can be 
ignored in our analysis. From effect (ii), we have that 



where nsDss dH is the SFD of asteroids with 2.5 < a < 2.6 AU contained in the 3rd. release of 
the SDSS-MOC, and n^now dH is the SFD of all the known asteroids in the same region, that can 
be computed from the catalog of proper elements contemporary to the SDSS-MOC. Figure [7] shows 
b{H) (full line) and its best fit (dashed line). The bias shows a linear dependence for H > 13.0, and 
a constant value ~ 0.22 for H < 13.0. 

The observed SFD ny"* is poorly known since, according to Table [TJ there are only two asteroids 
observed by the SDSS-MOC in the range 2.5 < a < 2.6 AU. Notwithstanding, we can use additional 
observations from the SDSS -MOC to get a more populated SFD. We follow the same approach as 



Roig and Gil-HuttonI (|2006l ) to identify V-type asteroids from the SDSS-MOC, but we disregard the 



information in the u band. This is justified since the u band is centered at ~ 0.35 ^m and it is not 
relevant for the identification of V-type asteroids. Thus, we considered the SDSS-MOC observations 
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with errors less than 10% only in the g, r, i, z bands and with any error in the u band. With this 
approach, we find 27 V-type candidates in the middle belt, including the three asteroids listed in 
Table [H 

The cumulative SFD of the 10 V-type candidates identified in the region 2.5 < a < 2.6 AU, 
and the corresponding debiased cumulative SFD computed from Eqs. ([6]) and are shown in 
Fig. [8] (full lines). It is worth recalling that both SFDs are affected by the completeness bias of 
the known asteroidal population with a > 2.5 AU. Therefore, we may compare these distributions 
to our predictions of Na>2.b computed using the observed Vesta family SFD which are affected 
by the same completeness bias. The comparison to Na>2.5 for Tage = 3.5 Gy (dashed line in Fig. 
[H|) indicates that the predictions are > 10 times smaller than required in the studied size range. 
However, we must bear in mind that the cumulative SFDs shown in Fig. [8] correspond to asteroids 
that are candidate V-type according to the SDSS colors. An unknown fraction of these bodies might 
not be true V-types, but belong to other taxonomic classes like 0-, Q-, R- or S-type. Therefore, 
our predictions could actually account for > 10% of the total population of V-type asteroids in the 
middle belt, assuming Tage ^ 3.5 Gy. 



6. Conclusions 

In this paper, we presented spectroscopic observations in the visible that confirm the existence 
of two V-type asteroids in the middle belt: (21238) 1995WV7 and (40521) 1999RL95. We investigate 
whether these two asteroids might have evolved from the Vesta family by slowly drifting in semi- 
major axis due to the Yarkovsky effect and crossing over the J3/1 mean motion resonance with 
Jupiter. Our results show that, in spite of the remarkable instability of the J3/1 resonance, km- 
size asteroids can cross it. The resonance crossing mechanism is probably not sufficiently efficient 
to explain the presence of all km-size V-type asteroids in the middle belt, but only some fraction 
(> 10%) of them. Most of these bodies either follow other p aths from the Vesta family or they 



come from a totally different source (e.g. ICarruba et al.l 120071 ). Notwithstanding, we cannot rule 



out the possibility that (21238) were a rare exception: the only one 5-km V-type asteroid in the 
middle belt that reached its present orbit by a very improbable, but not impossible, crossing over 
the J3/1 resonance. Only the discovery of more V-type asteroids in the middle belt, and the better 
knowledge of the SFD of these bodies may shed some light on this problem. 



This work has been supported through several grants and fellowships by the Brazilian Council 
of Research (CNPq), the NASA's Planetary Geology & Geophysics Program, and the Rio de Janeiro 
State Science Foundation (FAPERJ). 



- 14 - 



REFERENCES 

Asphaug, E.: 1997, 'Impact origin of the Vesta family'. Meteor. Planet. Sci. 32, 965-980. 

Binzel, R. P., G. Masi, and S. Foglia: 2006, 'Prediction and Confirmation of V-type Asteroids 
Beyond 2.5 AU Based on SDSS Colors'. ^4^45 / Division for Planetary Sciences Meeting 
Abstracts 38, 71.06. 

Bogard, D. D. and D. H. Garrison: 2003, '39Ar-40Ar ages of eucrites and thermal history of asteroid 
4 Vesta'. Meteor. Planet. Sci. 38, 669-710. 

Bottke, Jr., W. F., D. Vokrouhlicky, D. P. Rubincam, and D. Nesvorny: 2006, 'The Yarkovsky 
and Yorp Effects: Implications for Asteroid Dynamics'. Ann. Rev. Earth Planet. Sci. 34, 
157-191. 

Bus, S. J. and R. P. Binzel: 2002a, 'Phase II of the Small Main-Belt Asteroid Spectroscopic Survey. 
A Feature-Based Taxonomy'. Icarus 158, 146-177. 

Bus, S. J. and R. P. Binzel: 2002b, 'Phase II of the Small Main-Belt Asteroid Spectroscopic Survey. 
The Observations'. Icarus 158, 106-145. 

Carruba, V., T. A. Michtchenko, and D. Lazzaro: 2007, 'On the V-type asteroids outside the 
Vesta family. II. Is (21238) 1995 WV7 a former member of the Eunomia family?'. Astron. 
Astrophys., in press. 

Carruba, V., T. A. Michtchenko, F. Roig, S. Ferraz-Mello, and D. Nesvorny: 2005, 'On the V- 
type asteroids outside the Vesta family. I. Interplay of nonlinear secular resonances and the 
Yarkovsky effect: the cases of 956 Elisa and 809 Lundia'. Astron. Astrophys. 441, 819-829. 

Duffard, R. and F. Roig: 2007, 'Two new basaltic asteroids in the Outer Main Belt'. ArXiv e-print, 
astro-ph/0704.0230. 

Durda, D. D., W. F. Bottke, D. Nesvorny, B. L. Enke, W. J. Merline, E. Asphaug, and D. C. 
Richardson: 2007, 'Size frequency distributions of fragments from SPH/N-body simulations 
of asteroid impacts: Comparison with observed asteroid families'. Icarus 186, 498-516. 

Gladman, B. J., F. Migliorini, A. Morbidelli, V. Zappala, P. Michel, A. Cellino, C. Froeschle, H. F. 
Levison, M. Bailey, and M. Duncan: 1997, 'Dynamical lifetimes of objects injected into 
asteroid belt resonances'. Science 277, 197-201. 

Hammergren, M., G. Gyuk, and A. Puckett: 2006, '(21238) 1995 WV7: A New Basaltic Asteroid 
Outside the 3:1 Mean Motion Resonance'. ArXiv e-print, astro-ph/ 0609420. 

Harris, A. W.: 1979, 'Asteroid rotation rates II. A theory for the collisional evolution of rotation 
rates'. Icarus 40, 145-153. 



- 15 - 



Ivezic, Z., S. Tabachnik, R. Rafikov, R. H. Lupton, T. Quinn, M. Hammergren, L. Eyer, J. Chu, 
J. C. Armstrong, X. Fan, K. Finlator, T. R. Geballe, J. E. Gunn, G. S. Hennessy, G. R. 
Knapp, S. K. Leggett, J. A. Munn, J. R. Pier, C. M. Rockosi, D. P. Schneider, M. A. Strauss, 
B. Yanny, J. Brinkmann, I. Csabai, R. B. Hindsley, S. Kent, D. Q. Lamb, B. Margon, T. A. 
McKay, J. A. Smith, P. Waddel, D. G. York, and the SDSS Cohaboration: 2001, 'Solar 
System Objects Observed in the Sloan Digital Sky Survey Commissioning Data'. Astron. J. 
122, 2749-2784. 

Landolt, A. U.: 1992, 'UBVRI photometric standard stars in the magnitude range 11.5-16.0 around 
the celestial equator'. Astron. J. 104, 340-371. 

Lazzaro, D., C. A. Angeli, J. M. Carvano, T. Mothe-Diniz, R. Duffard, and M. Florczak: 2004, 
'S'^OS^: the visible spectroscopic survey of 820 asteroids'. Icarus 172, 179-220. 

Lazzaro, D., T. Michtchenko, J. M. Carvano, R. P. Binzel, S. J. Bus, T. H. Burbine, T. Mothc-Diniz, 
M. Florczak, C. A. Angeli, and A. W. Harris: 2000, 'Discovery of a Basaltic Asteroid in the 
Outer Main Belt'. Science 288, 2033-2035. 

Levison, H. F. and M. J. Duncan: 1994, 'The long-term dynamical behavior of short-period comets'. 
Icarus 108, 18-36. 

Marzari, F., A. Cellino, D. R. Davis, P. Farinella, V. Zappala, and V. Vanzani: 1996, 'Origin and 
evolution of the Vesta asteroid family.'. Astron. Astrophys. 316, 248-262. 

Masi, G., S. Foglia, and R. P. Binzel: 2006, 'Asteroid taxonomy in the Small Main Belt Asteroid 
Spectroscopic Survey system for 40313 minor planets fromthe 3rd Release of the Sloan Digital 
Sky Survey Moving Object Catalog'. Icarus, Submitted. 

McCord, T. B., J. B. Adams, and T. V. Johnson: 1970, 'Asteroid Vesta: Spectral Reflectivity and 
Compositional Implications'. Science 168, 1445-1447. 

Michtchenko, T. A., D. Lazzaro, S. Ferraz-Mello, and F. Roig: 2002, 'Origin of the Basaltic Asteroid 
1459 Magnya: A Dynamical and Mineralogical Study of the Outer Main Belt'. Icarus 158, 
343 359. 

Migliorini, F., A. MorbidelH, V. Zappala, B. J. Gladman, M. E. Bailey, and A. Cellino: 1997, 'Vesta 
fragments from v6 and 3:1 resonances: Implications for V-type NEAs and HED meteorites'. 
Meteor. Planet. Sci. 32, 903-916. 

MorbidelH, A., D. Nesvorny, W. F. Bottke, P. Michel, D. Vokrouhlicky, and P. Tanga: 2003, 'The 
shallow magnitude distribution of asteroid families'. Icarus 162, 328-336. 

Morbidelli, A., V. Zappala, M. Moons, A. Cellino, and R. Gonczi: 1995, 'Asteroid families close to 
mean motion resonances: dynamical effects and physical implications'. Icarus 118, 132-. 



16 



Moskovitz, N. A., M. Willman, S. J. Lawrence, R. Jedicke, D. Nesvorny, and E. J. Gaidos: 2007, 'A 
survey of basaltic asteroids in the Main Belt'. In: Lunar and Planetary Institute Conference 
Abstracts, Vol. 38. pp. 1663-1664. 

Mothe-Diniz, T., F. Roig, and J. M. Carvano: 2005, 'Reanalysis of asteroid families structure 
through visible spectroscopy'. Icarus 174, 54-80. 

Nathues, A., S. Mottola, M. Kaasalainen, and G. Neukum: 2005, 'Spectral study of the Eunomia 
asteroid family'. Icarus 175, 452-463. 

Petit, J.-M. and P. Farinella: 1993, 'Modelling the outcomes of high- velocity impacts between small 
solar system bodies'. Cel. Mech. Dyn. Astr. 57, 1-2. 

Roig, F. and R. Gil-Hutton: 2006, 'Selecting candidate V-type asteroids from the analysis of the 
Sloan Digital Sky Survey colors'. Icarus 183, 411-419. 

Rubincam, D. P.: 1995, 'Asteroid orbit evolution due to thermal drag'. J. Geophys. Res. 100, 
1585-1594. 

Rubincam, D. P.: 2000, 'Radiative Spin-up and Spin-down of Small Asteroids'. Icarus 148, 2-11. 

Tanga, P., A. Cellino, P. Michel, V. Zappala, P. Paolicchi, and A. dell'Oro: 1999, 'On the Size 
Distribution of Asteroid Families: The Role of Geometry'. Icarus 141, 65-78. 

Tedesco, E. F.: 1989, 'Asteroid magnitudes, UBV colors, and IRAS albedos and diameters'. In: 
R. P. Binzel, T. Gehrels, and M. S. Matthews (eds.): Asteroids II. pp. 1090-1138. 

Tholen, D. J.: 1989, 'Asteroid taxonomic classifications'. In: R. P. Binzel, T. Gehrels, and M. S. 
Matthews (eds.): Asteroids 11 pp. 1139-1150. 

Thomas, P. C., R. P. Binzel, M. J. Gaflfey, A. D. Storrs, E. N. Wells, and B. H. Zellner: 1997, 'Impact 
excavation on asteroid 4 Vesta: Hubble Space Telescope results'. Science 277, 1492-1495. 

Vokrouhlicky, D.: 1999, 'A complete linear model for the Yarkovsky thermal force on spherical 
asteroid fragments'. Astron. Astrophys. 344, 362-366. 

Vokrouhhcky, D., A. Milani, and S. R. Chesley: 2000, 'Yarkovsky Effect on Small Near-Earth 
Asteroids: Mathematical Formulation and Examples'. Icarus 148, 118-138. 

Wisdom, J.: 1982, 'The origin of the Kirkwood gaps - A mapping for asteroidal motion near the 
3/1 commensurability'. Astron. J. 87, 577-593. 

Zappala, V., A. Cellino, P. Farinella, and Z. Knezevic: 1990, 'Asteroid families. I - Identification by 
hierarchical clustering and reliability assessment'. Astron. J. 100, 2030-2046. 



This preprint was prepared with the AAS IfT^jX macros v5.2. 



-17- 



Fig. 1.— Visible spectra of (21238) and (40521) observed with CMOS at Gemini South (black 
lines). The gray lines are the spectra of other known V-type asteroids. The spectra are normalized 
to 1 at 5500 A and shifted by 0.5 in reflectance for clarity. A running box of 20 A has been applied 
to smooth the spectra. 
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Fig. 2. — Evolution of the semi-major axis, eccentricity and inclination of two test particles crossing 
over the J3/1 MMR at two different Yarkovsky drift rates. Note the different scales in the time 
axis. 
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Fig. 3. — Detail of the trajectory during the resonance crossing of the slowly drifting test particle 
shown in Fig. [2l The top panel shows the resonant angle a as a function of time. The bottom 
panel show the trajectory in the space a, a. The arrows indicate the direction of the trajectory. The 
numbers provide the temporal sequence of the trajectory. 
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Fig. 4. — The fraction /reach of V'esta's fragments that reached the J3/1 MMR border along three 
different values of Tage, and for two different values of Vey. 150 m/s (full lines) and 300 m/s (dashed 
lines). 
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Fig. 5. — The cumulative size distribution of the Vesta family after subtracting the known interlopers 
(full line) and three power law fits for different size ranges (dashed lines). Diameters were computed 
assuming albedo 0.4 for all bodies. 
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Fig. 6. — The cumulative distribution A^a>2.5 for three different values of r^ge and two models of 
the SFD: the observed SFD (full lines) and the single power law SFD (dashed lines). In all cases 
Vej = 200 m/s. The stars represent the known V-type asteroids in the middle belt. 
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Fig. 7. — The bias function (full line) in the region 2.5 < a < 2.6 AU, and its best fit (dashed 
line). The distributions nknow^-f^ (gray histogram) and Us^ss^H (outlined histogram) are shown 
for reference. Both distributions are normalized such that ^nAiJ = 1. 
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Fig. 8. — The cumulative SFD of V-type asteroids with 2.5 < a < 2.6 AU observed by the SDSS (full 
black line) and the corresponding debiased distribution (full gray line). The dashed line represents 
our predicted Na>2.5 for r^ge = 3.5 Gy, computed from the observed Vesta family SFD. 
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